Amplitude death state for hearing 
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We propose amplitude death phenomenon as an underlying mechanism of auditory transduction. 
When non-identical auditory hair bundles are elastically coupled, their spontaneous oscillations can 
be quenched to form an amplitude death state. We show, in this state, the hair cells are quiet 
and ready to detect oscillatory stimulus with coupling-strength dependent amplification. Numerical 
demonstration of the mechanism suggests that the non-uniformity of coupled hair cells can contribute 
to noise-robust auditory transduction. 



One of the surprising aspects of hearing is that the 
sound we hear may cause similar or smaller displacement 
of hearing organ than we expect from thermal fluctuation 
in equilibrium In ear, hair cells transform mechan- 
ical stimuli into neuronal signal with great sensitivity. 
Gold^2] has pointed out the ear has high damping due to 
the fluid in it and argued that the ear must contain active 
amplifier to achieve the high sensitivity and frequency se- 
lectivity. The hair cells indeed aniplify their inputs rely- 
ing on active feedback process [3, Gold's regenerative 
mechanism for hearing naturally lead to predict emitting 
sound of the ear due to the positive feedback oscillations. 
The ear indeed emits sound which is known as sponta- 
neous otoacoustic emission 5] and individual hair cells 
also show spontaneous motion of hair bundles The 
signal amplification might rely on the entrainment of the 
external signal to its spontaneous motion as usual radio 
engineering. However, direct application of the entrain- 
ment mechanism to hearing is problematic, because, if 
the hair bundles move spontaneously, the auditory neu- 
ron must receive strong noisy signal even in the absence 
of the external sound. 

The hair cell must be able to amplify weak stimuli and 
must be quiet in the absence of sound signal, which is 
coined here as quiet amplifier. The hair cells described 
by a generic mathematical modelfs^ and a biophysical 
model[9| can be quiet amplifiers. In these models, the 
hair cells are assumed to operate at the critical point of 
Hopf bifurcation [lol|, so the active amplification of weak 
stimuli is possible without relying on the entrainment 
mechanism. In reality, however, there exists noise in the 
process of gene expression thus it is more natural to as- 
sume a distribution of physical parameters of hair bun- 
dles rather than assuming all hair bundles are equivalent 
and located at the critical point. Furthermore, the ability 
to amplify weak signal is limited by temporal fluctuations 
which hides criticality. 

In this Letter, we provide a mechanism how the hair 
cells can be quiet amplifiers in spite of the noise and 
non-uniformity. We will show coupled hair bundles can 
enhance amplification and suppress mechanical fluctua- 
tions through amplitude death phenomenon which means 
quenching of the oscillation due to the coupling of the 
oscillator. This intriguing phenomenon was noted in the 



FIG. 1: The approximate phase diagram of two coupled oscil- 
lators when (sl) Hi — fi2 — i-L > and (b) /ii = — /i2 — fi > 0. 
In the regime of LOCKING, the natural frequencies of the 
two oscillators are locked with arbitrary relative phase, while 
the oscillation frequencies are different in the INCOHERENT 
regime. The amplitude of oscillation is zero in the amplitude 
DEATH regime. 



19th century by Rayleigh[ll[, who found that adjacent 
organ pipes can suppress each other's sound. Amplitude 
death is universal behavior that appears when any two 
or many different oscillators are mutually couDled |12l4l4| . 
The phenomenon can help auditory transduction as we 
will demonstrate in this work using numerical simulation 
of coupled hair bundles. While noisy fluctuations and 
unwanted spontaneous oscillation is strongly suppressed, 
the oscillators in the amplitude death state react sensi- 
tively to weak oscillatory stimulus. In this viewpoint, the 
disorder in hair bundles' size can contribute to hearing. 

A simple generic mathematical equations describing 
amplitude death is 



ii = (/ii + iwi - |zip)zi fc(z2 - Zl) 

Z2 = (M2 + «W2 - |Z2|^)^2 + - Z2), 



(1) 



where 2: is a complex variable of time, Wi is natural fre- 
quency of individual oscillator, and k represents the cou- 
pling strength. In the absence of the coupling {k = 0), 
the model describes weakly nonlinear system near a su- 
percritical Hopf bifurcation [l,[lH. For a single oscillator, 
the oscillatory solution Zi — y^e*'^'* becomes stable and 
the fixed point z = becomes unstable when fj,i > 0. 
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FIG. 2: The positions of three coupled hair bun- 
dles as a function of time. To investigate the effect 
of various coupling strength, the coupling constant k{t) 
changes as a function of time as shown in the inset ( 
in units of pN/nm). We choosed (/max[pN],fcsp[pN/nm]) = 
(337,0.58), (329,0.69), (357,0.58). Offset values for the dis- 
placement are added for clearance. 



When two oscillators have stable oscillatory solutions 
with different frequencies, coupling between the oscilla- 
tors can induce instability of the oscillatory motion. For 
instance, when = /ii = /i2 > 0, the amplitude death 
appears 1 13| if 



Auj >2n, and {I + { > - > 1, 

2 2/i /i 



(2) 



where Auj = wi — W2 is the frequency mismatch (Fig. 1 
(a)). Also if one of the oscillators has stable oscillatory 
solution, while the other is at its stable fixed point, the 
oscillatory motion can be quenched by the coupling. For 
instance, when ^ = fii — —fi2 > (Fig. 1(b)), the sta- 
bility of the fixed point zi = Z2 = is determined by the 
sign of real part of the Lyapunov exponent A satisfying 
the eigenvalue equation; ((^ -I- fc -I- A)^ -I- ujf){{^ — k — 
X)^ +iol) + 2k^{fi + k + X){p-k-X) + 2k'^LOiLU2 + k^ = 0. 
The amplitude death regime arises when both of the two 
eigenvalues have negative sign. This condition is met 
when 



|AiLj| > 2/i, and k > 



\Aio\fj. 



v/|Aw|2 - V 



(3) 



Note that in any cases the amplitude death requires cer- 
tain finite frequency mismatch Auj which was not usually 
assumed in previous hearing research. 

To investigate the response of the amplitude death 
state to oscillatory stimulus, let us consider the cou- 
pled oscillators subjected to periodic forcing as Zi — 
{fii + iuji — \zi\'^)zi + k{zj — Zi) + i^e'"*. The external pe- 
riodic force Fe^'^* excites the oscillators in the amplitude 
death state, and the oscillation frequencies are locked to 



the driving frequency oj. Analytic expression for the os- 
cillation amplitude and phase can be obtained when we 
consider a special case /i = /ii = /12 and the locking 
frequency w matches with the mean oscillator frequency 
Lj = {uji + ijj2)l2. Due to the symmetry, we can write 
z\(t) = zq cxp(iw<) and z^it) = Zq exp(ia;i). The com- 
plex amplitude zq satisfies a simple equation 

(M + i^-|^o|^)^o + fc(z*-zo)+i^ = 0. (4) 

The phase difference between the oscillators depends 
only on the amplitude of the oscillation and the coupling 
strength; 



-pi 



Aw 

2k = cot 0, 



(5) 



where zq = Rexjp{i(j)). It is interesting to note that there 
is intrinsic phase difference between the oscillators even 
for very weak stimulus limit. 

For weak stimulus, the phase difference is given by 
cos 6 ~ -p+2fc ^j^g amplitude is approx- 

imated as a linear response F cos 4>/{k{l — cos 2(/)) — 
fi). Nonlinearity becomes significant for i? > 1, where 
the amplitude R approximately satisfies R^ — fiR w F 
and the oscillators move in same phase cos0 w 0. In 
this case, the 1/3 law appears as in the critical oscilla- 
tor model 0. So, the mechanical response at the average 
frequency w = cj in the amplitude death state is approx- 
imated as 



R 
R 



1 

-F 



^(1_2|)2 + (AH)2 



(i?<l) 



(i? > 1) 



(6) 
(7) 



The oscillation amplitude R is linear in F for weak 
stimulus limit in contrast to the essential nonlinearity 
of a single oscillator of = ojH. As can be seen in 
Eq.([S]) and Eq.Q, the sensitivity defined hy R/F can be 
enhanced strongly as the coupling strength k increases, 
and even diverges at the verge of the amplitude death 
regime k — §(1 + (■^)^- Thus, the essential nonlinearity 
appears again at this special coupling. 

Now we investigate the possibility of the amplitude 
death phenomenon in a concrete biological system. For 
this purpose, we choose hair bundles of a bullfrog's hair 
cell, which has been concretely modeled in Ref. [^jli?]. Ac- 
tive hair-bundle movements results from Ca^~'"-dependent 
activity of the molecular motors which are connected 
to transduction ion channels through gate springs. We 
assume elastic coupling between hair bundles then the 
equations of motion for the position of i-th hair-bundle 
Xi and its molecular motor Xa.i are written as 



A— 

dt 



dt 



Fgs,i — ksp.iXi + k{Xi^i + ATi+i — 2Xi) 

-Fss,i + Fa,^, {i^l,2,---,N) (8) 
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where A and Aa are parameters relating the force and ve- 
locity. We used closed boundary condition Xq = X^+i = 

0, Xafi = ~ 0, FgsS ~ —kgs{Xi — Xa^i — DPo^i) 

is the force generated by the gating spring, Fa^i = 
0.14/max,i(l — 0.65Po,j) is the active force from the molec- 
ular motor, and Po,i — l/(l-l-exp( '~ /cVj^jf ' "'" )) is the 
open probability of the i-th ion channelQ. 

The non-uniformity of the hair bundles has been sim- 
ulated by introducing a distribution of the maximal gat- 
ing spring force /max and the bundles' pivot stiffness fcgp. 
Fig. 2 shows the positions of three hair bundles which are 
coupled with slowly time- varying constant k{t) as shown 
in the inset. One can see that the oscillation of three 
hair bundles becomes quenched simultaneously. Since 
the coupling constant for the quenching depends on de- 
gree of mismatch between hair bundles, we can not say in 
this work universal value of k for the bifurcation, but the 
appearance of the amplitude death is clear. The differ- 
ence between the k values for the quenching and revival 
of the oscillation indicates that the bifurcation is subcrit- 
ical. 

We will show the amplitude death mechanism survives 
in a more realistic model including effects of inertia, 
noise, and large number of hair bundles. Hair bundles 
of human's outer hair cells are coupled through tectorial 
membrane which has finite mass and stiffness. To mimic 
the inertia of the membrane, we introduce the coordi- 
nate of the elastically coupled massive elements Si which 
are also coupled to hair bundles (See Fig. 3 (a)). The 
governing equations of motions for the N coupled hair 
bundles are 

N 

mS'i = -7„iSi + fc(Si+i - 2Si + Si_i) + ^ ./HBj'5i^2j-i, 

j=i 

AXj = -/HB,j - fcgs(Xj - - DP,;,) - /c,p,jXj (9) 

where i = 1,2, 2N is index for the mass, j = 1, 2, N 
is the index for hair bundles, 7™ is the frictional constant 
of the mass. The j-th hair bundles is assumed to be tight 
to the 2j-l th mass, S2j-i = ^j, and /hbj is the force 
exerted on the mass exerted by j-th hair bundle. 

The inclusion of inertia, spatial dimensionality, or de- 
tail shape of the coupling is not a important factor to 
formation of the amplitude death state. The key to the 
death state is the non-uniformity of the hair bundles. 
Using Gaussian random number generators, we simulate 
random distribution of parameters around /max — 342 
pN and kg-p = 0.65 pN/nm with the variance about 50pN 
and 0.3 pN/nm, respectively. Fig. 3 shows numerical 
results from a particular set of the parameters which 
contains 10 self-oscillating hair bundles and 40 static 
bundles in the absence of coupling (Fig. 3(b)). As the 
coupling strength increases, the mutual coherence of the 
bundles' oscillation is developed. We find that at a cou- 
pling strength which is not unreasonable value fc w 8 
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FIG. 3: (a) Schematic for the elastically coupled hair bundles 
with massive membrane, (b) The stroboscopic view (snap 
shots at every 10 ms for 1 second.) of 50 uncoupled hair 
bundles, (c) The same view of amplitude death state {k = 10 
nN//jm). We used m=2 jUg. 

pN/nm, the oscillation become quenched and a stable 
amplitude death state appears (Fig. 3 (c)). 

The neural response of the hair cells to stimulus is 
not directly proportional to the mechanical displacement 
of the hair bundles, so the sum of displacement of all 
the bundles does not have direct significant meaning in 
auditory transduction. The influx of Ca?^ triggers the 
release of the neural transmitter, and also affects the con- 
ductance of channel thereby changes the membrane 
voltage of the hair cell. Therefore, the open probability 
of the transduction ion channel Pq is more relevant to 
the neural information. Note that multiple of hair cells 
are involved in auditory transduction for a single tone. 
(The human cochlear has ~ 3,500 inner hair cells and ~ 
11,000 outer hair cells. Inner hair cells are innervated by 
auditory nerve.) In this respect, the averaged value of 
the open probability 

1 ^ 

is encoded in auditory nerve potential. The sound stim- 
ulus is delivered to hair bundles through fluid velocity of 
the limp liquid in cochlea and also through movement of 
basilar membrane. The effect is modeled by adding ex- 
ternal force F{t) to the exerted on the hair bundle. For 
the noisy environment, we include also noise force i^N(0 
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FIG. 4: (a) The signal force applied on hair bundles as a 
function of time in units of pN. The ion-channel open proba- 
bility averaged over 50 hair bundles for (b) incoherent (fc = 0), 
(c) locking {k — 2pN/nm), and (d) amplitude death dominant 
state {k = 10 nN//im) of 50 hair bundles. The noise strength 
was chosen Dn = 2.8 10"^^ N s. 



i.e., we add F{t)+FNj{t), in the right hand side of Eq.® 
where < FNj{t)FN.j' (s) >= ^DfqSit — s)5j_ji {Dn is the 
noise strength). 

Fig. 4 shows, in the absence of signal ( when F{t) — 
in Fig. 4 (a)), the hair bundles movement have three 
classes. Incoherent regime exists where the averaged open 
probability has random fluctuations due to the distribu- 
tion of different hair bundles (Fig. 4 (b)). Locking domi- 
nant: The incoherent regime evolves to coherent locking 
regime where the hair bundles develop mutual coherence 
as the coupling increases |15| (Fig. 4 (c)). In this regime, 
due to its strong spontaneous oscillation in P* , it is diffi- 
cult to detect when the signal arrives. This problem can 
be overcome when the hair bundles are in amplitude death 
state (Fig. 4 (d)). In this state, fluctuations are strongly 
suppressed because other hair bundles' spontaneous mo- 
tions are quenched. One can see the noisy fluctuation in 
the absence of coupling ( Fig. 4 (b)) is suppressed in the 
amplitude death regime (Fig. 4 (d)). The sound signal is 
amplified and properly encoded in Pq so that it is clearly 
visible when the signal is arrived in Fig. 4 (d). 

In summary, we have proposed a mechanism for audi- 
tory transduction which operates robust non-uniformity 
and noise. When coupled non-identical oscillators form 
amplitude death state, oscillatory stimulus can be de- 
tected with active amplification and precise time-arrival 
information. The oscillation amplitude driven by oscil- 
latory stimulus is linear, but its sensitivity can be con- 
trolled by the coupling strength which can even diverge 
at the boundary of the amplitude death regime. Beyond 
the linear regime of the stimulus strength F, conventional 
nonlinear regime R ^ F^^^ appears showing compressive 
nonlinearity. We demonstrate the amplitude death mech- 



anism through realistic biophysical models for elastically 
coupled bullfrog's hair bundles, which indicates the non- 
uniformity of the hair bundles can contribute to auditory 
transduction through amplitude death mechanism. The 
mechanism proposed here is not sensitive to spatial di- 
mension or number of oscillators and it is worth investi- 
gating the auditory transduction in non- vertebrates, for 
instance, insects of which neuron spikes are coupled to 
active antennal motion [ig!]. 

Dierkes et.al.|l5| have studied coupled hair bundles 
where the coupling caused synchronization but didn't 
show amplitude death. It seems that the bundles went 
directly from incoherent to locking regime. We think 
the reason is probably due to their non-uniformity ( fre- 
quency gradient ) was not sufficient to cause the ampli- 
tude death (See Fig. 2 (a) and Ref.fii]). 

I thank Jung- Wan Ryu, Sergej Flach, Boris Alt- 
shuler, Martin Sieber, Jan Wiersig, Taegeun Song, Daniel 
Robert, Tom Duke and Pascal Martin for useful discus- 
sion. This work was supported by Basic Science Research 
Program through the National Research Foundation of 
Korea (NRF) funded by the MEST (2011-0009557). The 
author thanks Marie Curie Fellowship (2010-2011) of EU 
through University of Bath, and the hospitality of Max 
Planck Institute for Physics of Complex Systems during 
the author's visit. 



[1] W. Bialek, Annu. Rev. Biophys. Biophys. Chem. 16, 455 
(1987). 

T. Gold, Proc. R. Soc. London B 135, 492 (1948). 
J. Howard and A. J. Hudspeth, Neuron 1, 189 (1988). 

A. J. Hudspeth, Curr. Opin. Neurobiol. 7, 480 (1997). 
D. T. Kemp, Arch. Oto.-Rhino.-Laryngol. 224, 370 
(1979). 

P. Martin, A. J. Hudspeth, Proc. Natl. Acad. Sci. U.S.A. 
96, 14306 (1999). 

P. Martin, A. D. Mehta, and A. J. Hudspeth, Proc. Natl. 
Acad. Sci. U.S.A. 97, 3183 (2000). 

V. M. Egufluz, M. Ospeck, Y. Choe, A. J. Hudpeth, and 
M. O. Magnasco, Phys. Rev. Lett. 84, 5232 (2000). 

B. Nadrowski, Pascal Martin, and Frank Jiilicher, Proc. 
Natl. Acad. Sci. U.S.A. 101, 12195 (2004). 
S. H. Strogatz, Nonlinear Dynamics and Chaos, Addison 
Wesley Publishing Company (1994). 
A. Pikovsky, M. Rosenblum, and J. Kurths, Synchro- 
nization: A universal concept m nonlinear sciences Cam- 
bridge University Press, Cambridge, (2001). 
D. G. Aronson, G. B. Ermentrout, and N. Kopell, Phys- 
ica D, 41, 403 (1990). 

J.-W. Ryu et.al., J. Korean Phys. Soc, 55, 395 (2009). 
P. C. Matthews and Seven H. Strogatz, Phys. Rev. Lett. 
65, 1701 (1990). 

K. Dierkes, B. Lindner, and F. Jiilicher, Proc. Natl. Acad. 
Sci. USA 105: 18669-18674 (2008). 

K.-H. Ahn, D. Robert, Y. D. Chung, J. Lee, K.-J. Kim, 
J. W. Ryu, ( in preparation ) 



[9 
[10 

[11 

[12 

[13 
[14 

[15 

[16: 



